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Consider an electrically neutral spin 1/2 particle of mass M and a magnetic moment [io 
(a Neutron) in interaction with an infinite straight wire carrying a current I and located 
along the z-axis. The magnetic field generated by the wire is given by (we use units with 
c=l = h) 

5 or {-V,x,Q) m 
*- z \ x 2 + y 2y/2 W 

where x, y are cartesian coordinates in the plane perpendicular to the wire. The Hamiltonian 
of the particle is given by 

+ (2) 

The motion along the z-axis is free and will be ignored in the following. Thus we get a 
two-dimensional problem with 

H = h + 2I *\X+l*)w- (3) 

The problem formulated above has a dynamical symmetry From the form of the 
Hamiltonian there is one obvious constant of motion 

h = xp y - yp x + —. (4) 

Also, the 'Runge-Lenz vector' defined by 

Ai= ^ +JTm\p^ J ^, (< = 1,2) (5) 
satisfies the following commutation rules 

[J 3 , Ai] = i e l3 Aj, [A t , H] = 0, [A h A,] = -i H J 3 { ™ M)2 - (6) 
Thus for negative energy (— E) eigenstates of H we may put 

we get 



{Jot, Jp\ — i tap-yJ^ (8) 

where a, = 1, 2, 3. Thus J a generate the algebra of S0(3) in the space of negative energy 
eigenstates of H. For the space of positive energy eigenstates of H we change —E to +E in 
above and get the commutation rules of the SO(2,l) algebra. 

The SO (3) dynamical symmetry leads to an accidental degeneracy in the bound state 
spectrum. One can show that (E < 0) 

and 

H=- {2 't M) \ (10) 
so that the bound state spectrum is given by 

(2/^M) 2 



E j+ L, M - 2M(i+ i )2 ( U ) 



where j = ^ |, • • • and m = —j, — j + 1, • • • , j; which exhibits the (2j + 1) fold degeneracy 
explicitly. 

It may be mentioned that even if the bound state spectrum is deduced on the basis of 
dynamical symmetry it is difficult to use it for the construction of wave functions in the 
coordinate representation which have an ackward form. As in the case of the Coulomb 
problem the momentum representation allows us to use the dynamical symmetry in an 
elegant manner. 

As an example we formulate and solve the integral equation for the Green's function 
G(p,p') in momentum space. The calculation follows Schwinger's || technique for the 
Coulomb Green's function except for one important modification. 

To go over to the momentum representation we put 

?=i\7 P (12) 

and use 



p 2vr 
so that 



Vp ^-ln\p-p'\ = 5{p-p') (13) 



= f / *(p0i^4^v- (14) 

Vp 27TJ I 

The momentum space Green's function G(p,p') for energy E satisfies 

(E - §j)G(p,P') - ^ / W'^-p' GQ",?) = S(P- P") (15) 

where Oi{% = 1, 2, 3) are the well known Pauli matrices and e 12 = 1 = —£21, = 0. Assuming 
that 

E = -Sl (16) 
is real and negative (p > 0) we use the coordinates (n , raj) where 

ra = 2 , rij = ^ (17) 

and 

n 2 + n 2 = l. (18) 

It is convenient to treat (ra , n) as a three dimensional vector n. Here and in the following 
three dimensional vectors will be written without the vector symbol. The scalar product of 
two 3- vectors a and b will be written as a.b, where a = (00, 01, a 2 ). 

The element of area (solid angle) on the surface the unit sphere of equation (15) is 

dn\dn2 



Now 



M=i 



(19) 



dn x dn 2 = d 2 p J I — ) (20) 



where the Jacobian J is given by 



where the matrix elements of M are 

M - — 2piPj 

Since for a matrix N 

In detN = tr InN 

we get 

In det(l — M) = In n 

where we use the property 

M 2 = M(l-n ). 



Thus 



which gives 



Further 



We define 



, ( n-no = (^)V*'). 



[n-n') 2 = {n Q -n! Q ) 2 + {fi-fS) 2 

= (l + n )(l + ^) ^ 2 



and 



m O) = -±U-\n) W^f)} + U{n>) (30) 

We get 

r{Q, n')-v J D(n, q")t{q", n')dtt" = s(n - n') (31) 

where 

„ = ^ (32) 
Po 

which looks covariant under 3-dimensional rotations due to the notational artifice. 

The introduction of the matrices U (n) although permitted would appear to be redundant. 
However, if we omit them the form of D(Q,Q') would not display in a convenient fashion 
the dependence on the 3-vectors n and n', although the left hand side of equation (33) 
should do so. We determine U(n) by having D(Q,Q') satisfy a differential equation where 
its rotational invariance is obvious. 

We note that 

aA(p-p') a A jp-p') 

^a - 2"(1 +wo)(l + n ) = — p — . (33) 

{n — n'YpQ \P — P \ 

Also 

(*^ { ^^ = **'<?-?) 

= 2% f0±n^\ 2 5{n _ n ,y (34) 



Thus 



(a A y P )(l + no)*(l + n'^U{n)D lJ-\n ! ) = -tp S(Q - Of). (35) 



Let 



U(n) = J^U(n) (36) 



then 



tri( n )_*P2_(£ A y)u(n)D(n, n') = 6(n - n'). (37) 

1 + no 

The right hand side of equation (37) suggests the use of polar coordinates, hence if the 
polar angles of (n ,n) are 9, <fi; 

pi = p tan I - I cos(p, 

p 2 =p tan (^\sin(j). (38) 



We them have 



^ A ^ = -^-^ (39) 

where e p and are unit vectors. 

Thus the operator acting on D(Q,fl') in equation (39) involves at most first order dif- 
ferential operators dg and d^. If it were invariant under rotations we can make D(Q,Q') 
invariant too. The unique candidate in that case is the operator a + ba.L where a and b are 
constants. We choose U (n) so that this holds true. 

Since 

(T.L = -iio.e^de + o^d^ - cot(0) + a.e p 8$) (40) 

the choice 

U(n) = f(9, 0)exp {^^-a.e^ (41) 

maintains the - icr.e^de term (apart from other addicional terms without derivatives). 

Working out the derivatives it is easy to see that one must choose g(9, 4>) = which gives 

= a.L - ill' 1 a.e^dgU) + % U' 1 W (42) 

where in the bracketed terms the derivatives act only on U. 
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If we impose the condition that the last two terms are independent of the cr-matrices we 

get 

U(n) = + % -sinQ a. e<p (43) 

and then 

U- x {n)-^—aA \7 p U(n) = 1 +a.L. (44) 



1 + n 



Hence 



(i +a.L)D(n,n') = 6(n,rt). (45) 

Our choice of U(n) gives 



^-m 1 ^) (46) 

which allows us to calculate D(Q,Q') explicitly by using equation (32). We get 



D(n,n') = ± ( 1 -™™ ,) . (47) 

2?r (n - n'^ 2 



Although we know D(Q, Q') in a closed form it is convenient to rewrite it in terms of the 
eigenfunctions of the operator 1 +a.L. The eigenfunctions of the operators J 2 , L 2 , S 2 and J3 
where S = | and J = L + S are known to be Yj^(9, 0) where ± correspond to j = £ ± ~ 
respectively f|. Also 

( 1 + <r.L)Y±(6, 0) = ±(j + ]:)Y±(e, 0) (48) 



so that 



D(Q,Q') = E t^t 0) l&V, ^) " *) $)} ■ ( 49 ) 

jm ^ ' 2 



Going back to equation (33) we get 

/ 



which for v > has poles at 



2/M M 1 

y = — — = j + t;- (51) 

Po 2 



This gives the bound state energies which are same as before. The unnormalized bound 
state wave functions are ( - 1+ra 1 - ) 3/ 2 U^Y^JJ), <ft). 

Next, one can do the sum indicated in equation (52) to write an integral representation 
for r(fi, f2') which in turn leads to one for the Green's function G(p,p'). 

This is done by using the identities: 

3 + 1 1*^^ = 1*7^+ ,-, ,wf, (52) 



and using 

i=/V _1 dp (53) 
n jo 

valid for n > 1. Although each form of the identity leads to a different looking integral 
representation they are related through integration by parts. We concentrate on the second 
form. 

We have 

r(Q, Q!) = 5(Q, Q') + v f\p-"A + - p u &-)dp (54) 

Jo 

where 

A± = Zy - * *jm(M) ^mV, </>')■ (55) 



Using the explicit forms of Y^,(0, 0) one can show that 

A + = j-(d p p + a.L)$, 

A- = i-(-ia.L + a p )$ (56) 

where 

$ = ((i_p)2 + p ( n _ n ')2)-^ (57) 
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Calculating the derivatives we get 



A + = —(1 - p a.n a.n')$ 3 , 
Air 

A" = — (p - a.n a.n')&. (58) 

47T 



Thus 



r(Q, Q') = S(Q, H') + ^£ d p( - P^™ <r.ra') - p"(<7-™ a.n' - p) ) • (59) 

Substituting this in equation (31) we get the integral representation for G(p,p'): 

S(p-p') v 1 1 1 

G(p,p ) = — — — + 



E-T 8ttM E-T E-T \p-p'\ 2 
| [2M{E - ^±^) - ia 3 (p Ap') + ||p-p'| 2 ) 
(/(-„)- J(iz) + 7(1 + „)- 7(1 - v)) 
+ip a A (p - p')(I{y) + I(-u) + 7(1 + v) + 7(1 - v))} 

where po and v are given by equations (20) and (34). Also, 

m = ?d P — — r , (60) 

Jo {(1 - p) 2 + p(n - n') 2 } 5 

T=— T' = — (61) 
2M' 2M 1 J 

and 

{n ~ n) =-M (E-T)(E-Vy (62) 
In conclusion we have obtained an integral representation for the Green's function for a 
Neutron in interaction with a linear current. The results were deduced for the case E < 0. 
By analytic continuation to E > one can identify the scattering amplitude if the distorted 
free propagators can be identified. In contrast to the Coulomb problem it is not easy to solve 
the scattering problem in the coordinate representation by separation of variables. Although 
a general momentum transfer dependence of (p — p')~ 2 is indicated many details need to be 
worked out. This will be done elsewhere. 
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